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Abstract. We study the charge profile of a Cso-FET (field effect transistor) as used in the experiments of 
Schon, Kloc and Batlogg. Using a tight-binding model, we calculate the charge profile treating the Coulomb 
interaction in a mean-field approximation. At low doping, the charge profile behaves similarly to the case 
of a continuous space-charge layer and becomes confined to a single interface layer for doping higher than 
~ 0.3 electron (or hole) per Ceo molecule. The morahedral disorder of the Ceo molecules smoothens the 
structure in the density of states. 
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and electrical properties of surfaces, interfaces, thin films, and low-dimensional structures - 74.70.Wz 
FuUerenes and related materials 



> 
m 
m 

o 



I 

o 
o 



Surface doping of high quality organic crystals has re- 
cently been achieved by Schon, Kloc and Batlogg 
^. This breakthrough has led to novel devices, new su- 
perconductors and measurements of transport properties 
over unprecedented wide ranges of carrier concentration 
and temperature. A tunable field effect transistor (FET) 
is used to inject carriers near to the surface of the organic 
crystal. In this paper we report on calculations of the re- 
sulting charge profile. We chose to investigate the case of 
a Ceo-crystal which has a special interest in view of the 
recent observations of superconductivity with a high Tc 
(~ 52 K) in hole doped samples. 

The FET used in the work of Schon et al. is illustrated 
in Fig. |l|(a). In the present work, we consider a Cgo crys- 
tal (fee lattice) with a [001] plane parallel to the gate. 
Undoped Cgo is a semiconductor with a 2 eV gap. How- 
ever, when an electric potential is applied between the gate 
electrode and the source/drain electrodes, either electrons 
or holes accumulate on the interface between the Ceo crys- 
tal and the gate dielectric. This leads to a doping of the 
Ceo interface planes which allows current to flow parallel 
to the interface, between the source and the drain. The 
resistance of this channel can be measured as a function 
of temperature and doping. The FET can be represented 
as a planar capacitance with equal and opposite charges 
on the metal gate and the Ceo planes near to the interface 
as shown in Fig. 0(b). We will calculate the charge profile 
as a function of the total charge induced at the interface 
(or equivalently on the metal gate) rather then in terms 
of the potential applied to the gate. We limit ourselves 



Gate 
dielectric 
(AI2O3) 



Gate 

electrode Source and drain 
electrodes 




(b) 



Positively (negatively) 
/ charged gate 



^ Electron 
(hole) 
doped 
layers 



e-mail: swehrli@itp.phys.ethz.ch 



Fig. 1. (a) Schematic picture of the FET as used in the ex- 
periments of Schon, Kloc and Batlogg. Picture taken from 
(b) Model of the electronic system in the FET. 



to the case of an ideal planar interface with no steps or 
imperfections. 

We begin by introducing a tight-binding model to de- 
scribe the electronic structure of the Ceo crystal. The ef- 
fects of Coulomb interactions are treated within a mean- 
field (or Hartree) approximation. The resulting density 
profile gives rise to highly structured density of states 
(DOS). However we show that the inclusion of morahe- 
dral (or orientational) disorder leads to a broadening and 
suppression of this structure in the DOS. 
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The Ceo molecule, which has icosahedral symmetry, is 
almost spherical. Therefore, in the solid, it orders natu- 
rally in the close packed fee structure. In the following 
we will assume the (hypothetical) unidirectional structure 
where all Ceo molecules have the same orientation Q . The 
space group for this structure is Fm3m and its primitive 
cell contains one molecule. The actual low-temperature 
structure with space group PaS is more complicated and 
has a unit cell with 4 molecules. However, in the present 
work, we are interested in the overall charge filling of the 
layers, for which details of the band structure are unimpor- 
tant. The lattice constant of the cubic cell (containing two 
molecules) is a = 14 A for pure Ceo (|]- We introduce here 
the unit vectors e^;, and which span the cubic cell. 
If one considers only a [001] layer, then the Ceo molecules 
form a 2D square lattice. In this case the primitive cell 
has a side length of 5 = ^ = lOA which is the distance 
between two neighboring molecules. The 2D primitive cell 
is spanned by the unit vectors ei = "Tf^^^ ~'~ ^^'^ 
e,y). 
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The kinetic energy is well described by a tight-binding 
Hamiltonian with nearest neighbor (n.n.) hopping 0,^: 



H- 



TB 



E 

{nja,n' j' a') 



(1) 



where ta'a(<5) are the hopping integrals depending on the 
relative position 6 of neighboring molecules. The label n 
denotes the different layers starting with the n — 1 layer 
closest to the gate. Within a layer n, the different sites 
are labelled by the index j. Finally, the index a denotes 
the different orbitals of the Ceo molecule. The brackets () 
indicate a summation over n.n. only. Summation over spin 
degrees of freedom is implicitly assumed. Hamiltonian (^ 
is invariant under translations parallel to the z — Q plane 
and can partly be diagonalized in the 2D k-space: 



(2) 



where |n — n'| < 1 and k is the two-dimensional wave- 
vector parallel to the planes. In the case of the conduction 
band, the structure of the hopping integrals ta'a(5) are 
given by the ti„ symmetry of the LUMO (Lowest Unoccu- 
pied Molecular Orbital) which is threefold degenerate Q . 
The corresponding wave-functions can be chosen such as 
to transform as x, y and z under the icosahedral symmetry 
group which reduces the number of independent hopping 
integrals to 4. They are given in Table 0. The matrix ele- 
ments Hn'a' ,nai}^) m (|^) Can bc calculated explicitly which 
yields for intraplanar processes (n' — n), 



Hnx,nxO^) = 2txx[cOs{kib) + COs(fc26)] 
Hny,nvO^) = 2tyy[cOs{kib) -1- COs(fc2 6)] 

:,„y(k) = 2txy[cos{kib) - cos(fc26)] 
- (^'■^ = 2t^^[cos(fci6) + cos(fc2fo)] 



ny 

-^nz.nz (k 
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Table 1. Hopping integrals (in meV) for hopping in the 
5 — (110) direction of the [e^^eye^] coordinate system. The 
sign of txy changes for hopping in the (llO) direction. Hopping 
integrals in other directions follow from rotations around the 
threefold (111) axes. Numerical values were taken from and 
adjusted by a factor 1.11 for a lattice constant a = 14.0 A of 
pure Ceo- 



and for interplanar processes (n' = n + I), 

Hn+-lx,nx 

Hn+-iy,ny{^) = 2tzz cos[{ki-k2)b/2]+2txx cos[{ki+k2)b/2], 
Hn+iz,nzO<-) = 2txxCos[{ki-k2)b/2]+2tyyCOs[{ki+k2)b/2], 
Hr^.ix,nz(k-) = i2txysm[{ki ~ ^2)6/2], (4) 
Hn+iy^riz(k-) = i2txysm[{ki + k2)b/2]. 

The band structure and DOS of a single layer (interplanar 
processes are turned off) are shown in Fig. g. As a compar- 
ison, the 3D DOS is shown as well. The logarithmic van- 
Hove singularities are clearly apparent in the DOS. More- 
over, one observes a symmetric DOS. This is explained by 
the substitution k k + [j, ^) which leads to H —H 
for the one-layer Hamiltonian (H). It implies also that all 
bands at (^,5^) cross at zero energy. Coupling to neigh- 
boring layers breaks this symmetry. The structure of the 
hopping integrals for the valence band are somewhat more 
complicated due to the fivefold degeneracy of the HOMO 
(Highest Occupied Molecular Orbital) with h„ symmetry. 
They are described in detail in Ref. Q and similar matrix 
elements as in (H) and (|^) were calculated with parameters 
taken from 

The Coulomb interaction can be included by adding 
the following term to the tight-binding Hamiltonian ([|), 

Hcoul = ^ C^Gatc( (5) 

+ 2 E ^ee{nj,n'j')nnjann'j'a', 
nja^n' j' a' 

where Unja is the number operator. The first term is the 
potential due to the positive charge on the gate. The sec- 
ond term is the electron-electron Coulomb repulsion. Here 
we follow Antropov et al. by using a screened Coulomb 
interaction for electrons on different sites and an on-site 
interaction Uq for electrons on the same site. 



Vee{nj,n'f) = 



Uo 



e\Rn 



if 
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R 



■nj — ^^^n'j' 



(6) 



with e being the dielectric constant of Ceo and R„j de- 
noting the positions of the sites. Note that this interaction 
does not distinguish between different orbitals. The dielec- 
tric constant e and the on-site interaction Uq are easily 
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Fig. 2. (a) 

coordinates 
are r(0,0), 
spins. Solid 



-0.1 0.1 

Energy (eV) 



Ceo conduction band of the 2D square lattice. The 
of the high-symmetry points in the Brillouin zone 
P(f ,0),Q(f , f ). (b) DOS per Ceo and for both 
line: 2 dimensions. Dashed line: 3 dimensions. 



accessible in the literature. We use the experimental esti- 
mation of e = 4.4 by Hebard et al. Q and the theoretical 
value of J7o = 1 eV by Antropov et al. . Somewhat dif- 
ferent values of e = 3.66 and Uq = 1.27 eV are proposed 
from theoretical calculations by Pederson et al. . 

The full Hamiltonian including the Coulomb term is 
a complicated many-body problem. We solve it within 
the mean-field approximation. From symmetry, solutions 
should be invariant under translations parallel to the planes. 
The mean-field charge distribution and the resulting po- 
tential are then independent of the in-plane site position 
J and the wave- functions can be decomposed as 



1 



(7) 



„j is the inplane component 
Na is the number of sites per layer. For a fixed k, 



where v in & band index and R° 
of R„j 

the perpendicular part tpv\i{na) of the wave-function is the 
solution of a Schrodinger equation with a partly discrete 
spectrum {E^\s}- In the mean-field approximation, every 
occupied electronic i^k state contributes \'ipv\^{na)\^ /Nu 
electrons to each orbital a in layer n. The sum over all 
these contributions leads, at zero temperature, to the elec- 
tron distribution 



p(n) 



]f E El'/'-k(r.a)|', (8) 



where p{n) is the total number of electrons per site in 
layer n. The factor 2 comes from the spin degrees of free- 
dom. Within the mean-field approximation, the Coulomb 



MF 
Coul 



Umf {n) = C/Gatc {n) + ^ p{n') Vee (nj, n'f) . (9) 

n'j' 



The sum in the second term as well as the potential cre- 
ated by the gate are divergent. However, these divergences 
cancel when the overall system is neutral. In order to make 
this more explicit, we add and subtract a hypothetical 
potential created by uniformly charged layers whose to- 
tal charge equals the actual charge of the Cgo layers. The 
mean-field potential can then be rewritten as a sum of two 
terms 



UMFin) = Ucs,p{n) + C/corr(")- 



(10) 



The first term, J7cap, is given by the superposition of a 
series of planar capacitance formed by the negative charge 
—ep{n') on each layer n' combined with the corresponding 
positive charge +ep{n') on the gate. Choosing ?7cap(l) = 
yields 



(n,n')-l], (11) 



where it has to be distinguished whether the given layer n 
lies inside or outside the capacitance with charges ±ep{n'). 
Using £ — 4.4 and & = 10 A leads to a prefactor Aire^ / e\/2h - 
2.9 eV. The second part of the mean-field potential, J/coir, 
includes the corrections arising from the discreteness of 
the charge distribution and is given by 

C/corr(n) - ^ p(n')l/corr(|ri-n'|) with (12) 

n' 

o' 

where 6^ is the surface of the two-dimensional unit cell 
and R„ lies in the plane n. In equation ( p^ ) we expect 
the largest corrections for n' — n since the difference 
of the potential created by a uniformly and a discretely 
charged plane is only sizeable at short distance. Vcorr 
cannot be expressed analytically since it involves an in- 
finite sum over a two-dimensional lattice. The result from 
numerical summation is given in Table | The function 
Vcorr(^'^ = l'^ ~ '^'D is rapidly decreasing and we will 
only retain the An — term in the following. The correc- 
tion term then becomes 



(13) 



C/corr("-) = -pin) ( 3.9— - Uo 



In fact, the second part of this correction overestimates 
the on-site interaction Uq- The contribution of this term 
is proportional to p(n) which reflects the fact that an elec- 
tron feels its own mean-field. This is an artifact of the 
mean-field approximation. In order to improve our model. 
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Table 2. The function Vcoi 
summation. 



r(Zin) as calculated by numerical 



we require that the on-site interaction should only be effec- 
tive when there is more than one electron on the molecule 
(p(n) > 1), leading to the correction term 



Coir(,"-j 



-3.9— -hC/omax[0,p(n) - 1]. (14) 
eh 




Fig. 3. Relative charging of the first layers as a function 
of total charge in the system (per area of the 2D unit cell). 
Both axis are in logarithmic units. Solid line: Conduction band. 
Dashed line: Valence band. 



The mean- field Hamiltonian was solved self-consistently 
at zero temperature by numerical means. The mean- field 
potential ( |l0| ) was used with the correction term ( |l^ ) . The 
resulting charge profile as a function of total charge is 
shown in Fig. y. We have done the calculation for both, 
conduction and valence band, which led to similar results. 
In the following we focus on the conduction band. It can be 
seen that more than 98% of the total charge is confined to 
the first layer for doping higher than ptot = 0.3. Further- 
more the confinement to the interface increases with the 
total charge. The mean distance of the charge distribution 
from the interface is defined as zq = X]r!>i "^d p{n)/ ptot, 
where d — 7 A is the width of one layer. For ptot < 0.1, 
this mean distance is found to follow a power law given 

o — 1/3 

by 2:o[A] « 3.8ptot with ptot in units of particles per 
area of the 2D unit cell. This power law behavior with 
respect to ptot is identical to the case of a standard space- 
charge layer in a continuous medium [l^ , which predicts 

zo[A] = 3.02 (/9tot"i^/e)-i/^ = 5.08 p^^J^ , where the ef- 
fective mass (perpendicular to the interface) of the bulk 
conduction band minima (X(1,0,0) points) is estimated 
to be mz = 0.92 me [|[|^. However, it is important to 
notice that our microscopic calculation gives a different 
prefactor. In addition, zq tends to saturate for ptot > 0.1, 
a regime that cannot be understood in the continuum 
model. In order to test the role of the correction term, we 
repeated the same calculation with C/mf = f^Cap, i-e. for 
uniformly charged planes. This led to a somewhat weaker 
confinement at high doping (~ 10% less charge on the first 
layer for ptot = !)■ At low doping the effect of C/coir van- 
ishes because the charge is distributed over several layers. 
Knowing the charge profile and hence UyiFin), the band 
structure can be calculated. Fig. ^ shows the conduction 
band for ptot = 0.1. Few bands lie below the continu- 
ous spectrum and only the lowest of the 3 LUMO-bands 
is occupied in the region of the Q-point. In Fig. |^(a) the 
evolution of the discrete states at Q is shown as a function 
of doping. One observes that the Fermi level is always be- 
low the bottom of the continuous spectrum. This is true in 
general and reflects the fact that all electrons are bound 
to the interface since the overall system is neutral. Fur- 
thermore, only the lowest band is occupied at low doping 
which is again consistent with theory of space-charge lay- 



0.05| 




Fig. 4. Solid lines: Conduction bands below the continuous 
spectrum for ptot = 0.1. Bold solid line: Onset of the continu- 
ous spectrum. Dashed line: Fermi energy. 
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Of course, at doping higher than ptot ~ 2 the 
second LUMO-band has to be filled as well. In addition, 
we note that the continuous spectrum shifts up sharply 
as the total charge increases. Therefore, at doping higher 
than Ptot ~ 0.3, the first layer becomes essentially decou- 
pled from the continuous states of the subsequent layers. 
The energy levels at Q remain constant and the dispersion 
is given by the LUMO-bands of a single, isolated layer 
(Fig. 0(a)). The decoupling effect is due to the correction 
term ([l4|) which drastically shifts the energy of the first 
layer with respect to the second one. If this term is omit- 
ted, as shown in Fig. ^(b), then the bottom of the contin- 
uous spectrum follows roughly the Fermi energy and the 
decoupling of the first layer is much less effective. 

In A3 Geo, a scaling of the superconducting Tc with 
increasing lattice constant and hence with the DOS is 
observed However, the sharp peaks of the 2D DOS 
(Fig. ^(b))are not reflected in the behavior of Tc as mea- 
sured by Schon et al. [Q. A possible explanation for this 
is the morahedral disorder which is considered in the fol- 
lowing for the case of a single 2D layer. 

The rotational motion of a Cgo molecule in a crystal 
is described by a potential in the three Euler angle co- 
ordinates. This potential has, in addition to an absolute 
minimum, a local minimum which is only A = 11 meV 
higher and corresponds to a non-equivalent orientation 
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Fig. 5. (a) Solid lines: Discrete levels at Q as a function 
of ptot- Bold solid line: Bottom of the continuous spectrum. 
Dashed line: Fermi energy, (b) As (a), but without correction 
term Ucon- 



of the molecule |Q. Therefore, Ceo molecules can flip 
into the second orientation causing the so-called mora- 
hedral disorder. The ratio between the number of flipped 
and unflipped molecules is given by a Boltzmann factor 
exp(— Zi/fcT). However, a freezing in is observed at 85 K. 
Below this temperature the crystal is in a glass phase with 
15% of the molecules being flipped. If the problem is con- 
sidered from a tight-binding point of view, then one ex- 
pects the hopping integrals to be different for hopping be- 
tween inequivalently oriented molecules. Hence, one gets a 
tight-binding model where different hopping integrals are 
distributed randomly. This will suppress structure in the 
DOS, as was calculated by Gelfand and Lu in the case 
of bulk Cgo |l^. Here, we consider the effect of morahe- 
dral disorder on the DOS of a single electron-doped [001] 
layer. The two types of inequivalent orientations lead to 
two sets of LUMO wave-functions. We shall denote the 
new ones on flipped molecules as \x), \y) and \z). These 
two sets can be related by a 90-degree rotation around 
the z-axes. The hopping integrals between inequivalently 
oriented molecules are given in Table In order to calcu- 
late the DOS of the disordered system we set up a finite, 
two dimensional system which could be diagonalized ex- 
actly (for non- interacting electrons). The resulting DOS 
is shown in Fig. ^. The van-Hove singularities are washed 
out, but a rather sharp peak remains at the center of the 
band. 

In conclusion, we have investigated the electronic prop- 
erties of surface doped Cgo. The charge profile was calcu- 
lated as a function of total charge in the system and it 
was found that the charge accumulation on the first layer 
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Table 3. Hopping integrals (in meV) for hopping from an 
unflipped to a flipped molecule in the (110) direction of a fee 
lattice. The new (tilded) wave-functions were ordered accord- 
ing to their parity. Note that the sign of t^x and tyy change for 
hopping in the (110) direction. Numerical values were taken 
from lis] and adjusted by a factor 11 as proposed by Gelfand 
and Lu. 




Fig. 6. DOS for a disordered 50x50 site lattice and a flip- 
ping ratio of 15%. The statistical average over 200 disorder 
configurations was taken. 



increases with the total charge. At high doping, above 
~ 0.3 electron per molecule, the first layer becomes es- 
sentially decoupled from the subsequent layers. This sug- 
gests that the electronic system should be well described 
by a single layer in this regime. In particular, the rele- 
vant DOS is then the one of a two dimensional system 
and hence substantially higher than in the bulk. Having 
chosen the unidirectional structure (one orientation per 
primitive cell), we found maximal values of the DOS near 
the center of the band. Introducing morahedral disorder 
led to a smoothened DOS which still shows a peak at half 
filling, i.e. in the region where Schon ct al. measure the 
maximum Tc- 
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